第 39 卷 第 4 期 天 x 学 4t JE Vol.39, No.4 
2021 年 12 月 PROGRESS IN ASTRONOMY Dec., 2021 
doi: 10.3969/j.issn.1000-8349.2021.04.02 
= ba SIN 进 
宇宙 黎明 和 再 电离 时 期 探测 进展 
EXEAT, APO (812, KARS, 3B RU, aE 

(1. 中 国 科 学 院 上 海天 文 台 ， 上 海 200030; 2. 中 国 科 学 院 射 电 天 文 重点 实验 室 ， 南 京 210033; ”3. 中 国 科 

学 院 大 学 ， 北 京 100049) 

摘要 : ?黎明 和 再 电离 信号 探测 是 目前 运行 以 及 即将 建成 的 低频 射电 阵列 (如 平方 公里 阵列 


square kilometer array, SKA) 
迁 产 生 的 21 cm 辐射 是 进行 宇 
刘 号 可 以 描绘 册 
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的 首要 科学 目标 2 


TS 
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黎明 和 再 电离 时 期 探测 的 最 
代 发 光 天 体 的 产生 和 形成 过 程 ， 从 而 研究 宇宙 


EE 


期 的 HI 超 精细 能 级 跃 


接 而 有 效 的 工 


L ， 通 过 探测 HI 的 


要 介绍 近年 来 宇宙 黎明 和 再 电离 探测 方向 的 相 


Xx 键 


中 图 分 类 号 : P159, P161 


探测 时 期 探测 中 关于 前 景 去 除 、 数 据 处 理 方法 等 相关 问 
量 、 功 率 谱 测量 方面 的 研究 进展 。 
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从 黑暗 走向 光明 的 
宇宙 黎明 和 再 电离 


ai 


镜 在 全 天 总 功率 测 


其 探测， 射电 干涉 阵列 


宇宙 大 爆炸 以 后 ， 当 宇宙 年 龄 为 38 万 年 时 ， 原 子 核 与 电子 复合 成 为 HI 原子 ， 光 子 得 


以 自由 地 运动 ， 从 而 形成 J 
发 光 天 体 形成 。 宇 宙 第 一 代 恒星 和 星系 开始 形成 ， 它 们 发 出 的 紫外 光 辐 


使 得 整个 宇宙 开始 明亮 起 来 ， 宇 宙 进 入 黎明 和 再 电离 时 期 。HI 21 cm 
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黎明 和 再 电离 时 期 的 21 cm 信号 红 移 到 今 


宇宙 背景 辐射 。 之 后 宇宙 进入 了 漫长 的 黑暗 时 期 ”"， 直 到 第 一 代 
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层 析 方 法 以 外 ， 男 一 种 方法 是 “21 cm 和 森林” 观测， 这 种 观测 是 以 高 红 移 ( 红 移 6 以 上 ) 的 
类 星体 或 伽 马 暴 余辉 等 一 些 非常 强 的 射电 源 作为 背景 辐射 源 ， 探 测 视 线 方 向 上 的 HI 在 它们 
各 自 的 光谱 上 产生 的 吸收 线 信号 ""。 不 同 红 移 的 结构 在 类 星体 光谱 的 不 同 频 率 处 产生 吸收 
线 ， 形 成 “森林 ” 似 的 光谱 结构 。 不 同 于 21 cm 层 析 方 法 的 是 ，21 cm 森林 信号 更 加 敏感 于 
星系 际 介质 的 温度 ， 能 够 更 有 效 地 提取 宇宙 温度 演化 的 信息 。 

正如 Gunn 和 Peterson" 所 预测 的 ， 在 两 个 z e 6 的 QSOs 中 检测 到 中 性 IGM (inter 
galactic matter， 星 系 际 物质 ) 对 Lya 的 吸收 ， 彻 底 改 变 了 对 宇宙 再 电离 时 期 的 理解 。21 cm 
森林 的 方法 ”可 以 提供 一 些小 规模 结构 的 信息 ， 例 如 作为 单个 暗物质 迷你 量 ， 小 质量 星系 ， 
甚至 早期 宇宙 的 了 H I 末 恒星 区 。 通 过 大 规模 的 数值 模拟 ，Carilli 等 人 ”发 现 一 些 高 密度 区 域 
Lya 森林 的 高 红 移 源 在 射电 品类 星体 光谱 中 产生 了 可 识别 的 吸收 特征 ， 但 是 他 们 无 法 解决 
最 小 质量 尺度 上 的 气体 云 替 塌 ， 尽 管 这 些 云 代表 了 高 红 移 IGM 中 最 丰富 的 团 块 。Furlanetto 
和 Abraham” 通 过 半 解 析 探 索 这 些 紧 凑 型 吸收 系统 的 统计 数据 ， 补 充 了 Carilli 等 人 的 研 
究 结果 。 这 种 半 解 析 方 法 能 够 检验 结果 对 输入 参数 变化 的 敏感 性 ， 可 以 更 好 地 通过 射电 吸收 
光谱 测量 物理 量 。Xu EA 以 高 红 移 点 源 作为 背景 ， 计 算 了 单个 迷你 暗 党 和 矮星 系 的 吸收 
信号 ， 包 括 内 部 和 周围 的 气体 含量 。 本 文 主要 讨论 层 析 法 的 进展 。 

然而 ， 宇 宙 黎 明和 再 电离 时 期 的 HI 21 cm 信号 非常 微弱 ， 相 对 于 宇宙 微波 背景 在 亮 温 
度 上 的 起 伏 的 最 大 值 只 有 约 0.01 K”， 比 包含 银河 系 辐射 以 及 河 外 射电 源 的 前 景 污染 低 4~5 
个 数量 级 ， 因 此 微弱 的 宇宙 黎明 和 再 电离 信号 被 强大 的 宇宙 前 景 信号 淹没 ， 要 将 前 景 污染 信 
号 移 除 到 5 个 数量 级 以 下 ， 在 目前 的 低频 射电 数据 处 理 中 还 存在 困难 。 另 外 ， 在 低频 射电 波 
段 存在 着 许多 人 为 和 干扰， 如 电视 基站 、 手 机 、 广 播 等 通信 干扰 ， 在 地 球 上 能 够 建立 低频 射电 
望远镜 的 射电 宁静 区 域 非常 稀少 。 总 之 ， 在 实际 探测 中 存在 着 诸多 困难 和 挑战 ， 例 如 前 景 的 
识别 和 去 除 ， 高 动态 大 视 场 成 像 ， 校 准 和 电离 层 改 正 ， 仪 器 效应 去 除 ， 海 量 数据 处 理 等 。 
目前 宇宙 黎明 和 再 电离 时 期 的 探测 在 研究 方法 上 主要 包括 三 个 方向 : 宇宙 黎明 和 再 电 
离 信号 的 全 天 总 功率 测量 ， 宇 宙 黎 明和 再 电离 信号 的 统计 测量 ， 即 功率 谱 测 量 ， 直 接 对 再 
电离 区 域 实施 成 像 观测 。 在 这 三 个 研究 方向 中 ， 全 天 总 功率 测量 通过 单 天 线 进行 观测 ， 宇 
宙 黎 明和 再 电离 信号 的 功率 谱 测 量 和 成 像 研究 则 是 基于 低频 射电 干涉 阵列 的 观测 来 实现 的 。 
宇宙 黎明 和 再 电离 信号 的 功率 谱 测 量 和 成 像 研究 是 目前 运行 的 和 即将 建成 的 低频 射电 阵列 
的 首要 科学 目标 之 一 。 通 过 不 同 红 移 处 的 HI 21 cm 辐射 观测 ， 可 揭示 不 同时 期 宇宙 中 的 HI 
= 分 布 ， 从 而 描绘 出 第 一 代 恒 星 和 星系 的 形成 及 气体 再 电离 的 复杂 过 程 ， 精 确 测量 宇宙 大 尺度 
吉 构 的 初期 演化 ， 为 解决 宇宙 结构 起 源 等 问题 提供 重要 信息 和 线索 。 
近年 来 ， 为 了 拨 开 迷雾 看 到 “宇宙 的 第 一 缕 曙 光 ”， 推 进 宇 宙 学 研究 进入 一 个 新 的 黄 
时 代 ， 字 宙 黎 明和 再 电离 时 期 探测 成 为 目前 国内 和 国际 上 的 低频 观测 设备 以 及 即将 开始 建 
设 的 大 型 射电 阵列 的 首要 科学 目标 之 一 。 进 行 全 天 总 功率 测量 的 实验 主要 包括 : EDGES 
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proposed dark ages radio explorer)” a , PRI“M (probing radio intensity at high-Z from 
marion) = , ASSASSIN (the all-sky signal short-spacing interferomete) = , DAPPER (dark 
ages polarimeter pathfinder)” 和 DSL (discovering the sky at the longest wavelengths with 
small satellite constellations) 等。 进行 宇宙 黎明 和 再 电离 时 期 统计 测量 的 低频 射电 阵列 主 
要 包括 : 21CMA (21 Centimeter Array)”, GMRT (Giant Meterwave Radio Telescope), 
MWA (Murchison Wide-field Array), LOFAR (Low Frequency Array) ^ ^, PAPER 
(Precision Array to Probe Epoch of Reionization) ™ ^! #1 HERA (Hydrogen Epoch of Reion- 
ization Array) 等， 这 些 目 前 已 经 建成 的 低频 射电 阵列 是 未 来 即将 建成 的 大 型 射电 阵列 
SKA 的 探 路 者 阵列 。 未 来 的 SKA 在 灵敏 度 和 分 辨 率 上 的 优势 将 具备 对 宇宙 黎明 和 再 电离 时 
期 进行 功率 谱 测量 和 直接 成 像 研究 的 能 力 。 本 文 对 近年 来 宇宙 黎明 和 再 电离 探测 领域 的 研 
究 进展 进行 介绍 。 

本 文 结构 如 下 : 第 2 章 概述 宇宙 黎明 和 再 电离 探测 在 观测 和 数据 处 理 过 程 中 的 相关 问 
Bi. 第 3 章 介 绍 目 前 低频 实验 的 探测 进展 ; 第 4 章 给 出 总 结 和 展望 。 


2 宇宙 黎明 和 再 电离 时 期 探测 概述 


2.1 21cm 信号 

HI 的 21 cm 辐射 是 由 HI 超 精细 能 级 跃迁 产生 的 。HI 的 电子 存在 着 自 旋 平行 和 反 和 平行 
两 种 状态 ， 当 这 两 种 状态 发 生 转 换 时 ， 便 发 生 了 超 精 细 能 级 间 的 跃迁 ， 产 生 了 波长 为 21 cm 
的 辐射 。 因 为 HI 发 生 自 发 超 精 细 能 级 跃迁 的 概率 非常 小 ， 对 于 1 个 五 原子 而 言 ， 发 生 一 次 
自发 跃迁 需要 的 时 间 是 大 约 1000 万 年 。 第 一 代 发 光 天 体 辐射 出 的 光子 加 热 HI， 使 其 周围 的 
HI 发 生 了 超 精 细 能 级 跃迁 。 

HI 的 21 cm 线 是 一 条 具有 特定 频率 的 谱 线 ， 在 不 同 的 探测 波段 对 应 于 不 同 红 移 时 期 的 
HI 信号 。 随 着 宇宙 的 膨胀 ， 宇 宙 歼 明和 再 电离 时 期 的 HI 21 cm 信号 传播 到 今天 的 波长 范围 
属于 米 波 波段 ， 也 就 是 低频 射电 波段 。 因 此 ， 可 以 在 低频 射电 波段 通过 探测 HI 21 cm 信号 
来 研究 第 一 代 发 光 天 体 的 形成 ， 从 而 研究 宇宙 的 黎明 和 再 电离 时 期 。HI 超 精 细 能 级 跃迁 产 
生 的 21 cm 信号 ， 可 以 通过 亮 温 度 来 描述 ， 表 示 为 ”: 


oTw(r, z) = To(z)[1 + bm(r)|[1 — zi(r)]. ， (1) 


fr 


RP, Om 是 物质 密度 ，z; 是 气体 的 电离 度 (o; = 1 表示 气体 完全 电离 ，z; = 0 表示 气体 是 
Sé PIER), To(z) 表示 为 : 


T; -T,N f Qh? 0.143\ /14- 2M]? 
To(z) — 2 > 
o(2) 6( Ts ) (5) 5 (G2) ( 10 ) l (2) 
HF, Ts æ HINA Weim, T, æ CMB 的 温度 ，h 是 普 朗 克 常 数 ，02 ETER, Rn 是 


总 物质 密度 ，z 是 红 移 。 假 设 Ts 在 整个 IGM 中 都 与 气体 温度 Tus 耦合 ， 并 且 和 气体 温度 远大 
T CMB 温度 。 一 旦 建立 了 充分 的 Lya 背景 以 将 自 旋 温度 与 气体 温度 耦合 ， 并 且 早 期 又 射 
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线 或 其 他 来 源 提 供 了 所 需 的 热量 ， 该 假设 便 成 立 。 当 HI 的 自 旋 温 度 高 于 5 
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的 亮 温度 时 ， 会 产生 21 cm 发 射 信号 ; 当 HI 原子 的 自 旋 温度 低 于 宇宙 微波 背景 辐射 的 亮 温 


度 时 ， 则 会 产生 21 cm 吸收 信号 。 因 此 可 以 通过 探测 相对 于 宇宙 微波 


起 伏 来 探测 HI 的 发 射 和 吸收 信和 号 。 
2.2 HERRAR 


在 利用 HI 21 cm 谱 线 重 构 宇 宙 演 化 史 的 过 程 中 ， 提 取 宇 宙 歼 明和 月 


题 之 一 就 在 于 能 否 准确 地 扣除 银河 系 和 银河 系 外 的 射电 前 景 污 染 。 前 景 主要 包括 射电 蛇 、 银 


Ab EL 
[za 


辐射 亮 温度 上 的 


了 电离 信号 的 关键 问 


河 系 同步 辐射 、 银 河 系 自由 -自由 辐射 、 银 河 系 外 射电 源 。 在 各 种 污染 性 前 景 成 分 中 ， 银 河 
弥漫 辐射 (包括 同步 辐射 和 自由 -自由 辐射 ) 和 银河 系 外 点 源 是 主要 污染 成 分 "一 ” ， 据 估计 ， 
近年 来 ， 国 际 上 在 前 


在 约 150 MHz 时 ， 它 们 可 能 分 别 约 占 前 景 污染 总 量 的 71% 和 2790". 


景 扣除 方面 开展 了 大 量 研究 ， 目 前 对 于 前 景 去 除 通 常 使 用 的 方法 基 


FA 


电 前景 谱 在 频率 空 
间 平 滑 性 的 假设 ， 通 过 使 用 光滑 函数 ”、 多 项 式 拟 合 ”、 主 成 分 分 析 (PCA)”、 奇 异 值 分 


解 (SVD)”、K-L 变换 ”等 方法 在 频率 空间 进行 射电 前 景 的 拟 合 和 去 除 。 然 而 ， 由 于 背景 
信号 比 前景 低 5 个 数量 级 (如 图 四 所 示 )， 因 此 将 前 景 通过 拟 合 精确 扣除 到 5 个 数量 级 以 下 还 
存在 着 困难 。 并 且 ， 由 于 仪器 系统 误差 的 存在 ， 在 观测 过 程 中 前 景 辐射 的 极 化 成 分 可 能 会 泄 
露 到 非 极 化 的 21 cm 信和 号， 增益 在 频 域 的 变化 等 都 会 影响 前 景 去 除 的 精度 ， 从 而 为 提取 宇 


宙 黎 明和 再 电离 信号 造成 困难 。 


I(-11)C(2x)?/10?K? 


1055 10 10° 


ik: 红色 细 实 线 为 点 源 ， 红 色 细 点 线 为 银河 系 同 步 辐射 ， 红 色 细 虚线 为 银河 系 外 


- 色 细 长 虚线 


为 银河 系 自 由 辐射 ，CMB 用 绿色 的 点 划 线 表示 。 


1 2=9.2 (v = 140 MHz) 时 的 21 cm 信号 功率 谱 ( 蓝 色 粗 线 ) 及 各 前 景 成 分 功率 谱 加 
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前 景 污染 中 ， 源 于 银河 系 的 污染 包括 银河 系 的 同步 辐射 和 自由 -自由 辐射 ， 比 宇宙 再 电 
离 背 景 信号 高 5 个 数量 级 。 由 于 在 频率 空间 认为 银河 系 的 辐射 谱 是 平滑 的 ， 因 此 可 以 通过 拟 
合 平滑 函数 的 方法 来 进行 扣除 。 河 外 射电 源 是 前 景 污染 的 主要 成 分 之 一 ， 其 中 形状 规则 的 点 
源 可 以 通过 高 斯 函数 拟 合 等 方法 进行 扣除 ， 但 是 具有 大 尺度 弥散 结构 的 展 源 却 是 前 景 去 除 
中 的 难点 。 前 景 污染 中 弥散 辐射 的 成 分 主要 包括 在 低频 射电 波段 具有 弥散 辐射 的 大 质量 星 
系 团 以 及 纤维 状 的 宇宙 大 尺度 机 构 ， 在 形态 上 与 宇宙 再 电离 背景 的 观测 尺度 类 似 ， 分 布 比较 
杂 ， 呈 非 高 斯 分 布 ， 所 以 很 难 通过 高 斯 函数 进行 精确 的 拟 合 。 根 据 MWA GLEAM™ “在 
低频 波段 70 ~ 230 MHz 的 巡天 结果 ， 探 测 到 了 200 余 个 射电 源 具 有 大 尺度 的 弥散 结构 。 随 
着 望远镜 性 能 的 提高 和 数据 处 理 方法 的 完善 ， 更 多 具有 弥散 辐射 的 射电 源 将 会 被 探测 到 。 这 
些 具 有 弥散 结构 的 射电 源 可 以 分 为 射电 晕 、 微 射电 晕 和 射电 遗迹 。 由 于 其 陡 谱 的 性 质 ， 在 低 
频 波段 呈现 出 更 强 的 辐射 ， 因 此 在 宇宙 黎明 和 再 电离 探测 的 频率 范围 内 ， 弥 散射 电源 的 影响 
是 不 能 被 忽视 的 。 尽 管 相 对 于 银河 系 前 景 ， 星 系 团 星系 际 介质 射电 展 源 的 辐射 较 弱 ， 但 仍 比 
宇宙 黎明 和 再 电离 时 期 的 HI 21 em 信号 强 2 ~ 3 个 数量 级 ”( 如 图 团 所 示 )。 
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银河 系 外 前 景 
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同步 辐射 
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N md oO prz >= [Eu 
图 2 120 MHz 低频 波段 的 宇宙 再 电离 21 cm 信号 及 各 类 前 景 示 意图 


目前 弥散 源 拟 合 方法 包括 多 高 斯 函数 拟 合 、shapelets 等 。shapelets 是 一 系列 由 加 权 
Hermite 多 项 式 组 成 的 正 交 基 函 数 ”。shapelets 已 在 许多 天 文 应 用 中 使 用 ， 包 括 模拟 尘埃 
的 三 维 分 布 ”、 在 射电 图 像 模拟 中 的 弱 引 力 透镜 测量 ”、 引 力 透 镜 图 像 ” 等 。 在 宇宙 黎明 
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和 再 电离 探测 中 ， 不 仅 需 要 在 像 空 间 进 行 前 景 成 分 的 拟 合 和 去 除 ， 还 需要 在 UV 空间 对 非 格 
点 化 的 数据 进行 操作 ， 直 接 进 行 功率 谱 的 提取 。 因 此 ， 目 前 弥散 源 的 拟 合 仍 然 是 宇宙 黎明 和 
再 电离 探测 的 难点 。 
2.3 ”数据 处 理 

由 于 宇宙 黎明 和 再 电离 信号 非常 微弱 ， 因 此 在 观测 和 数据 处 理 的 各 个 环节 都 有 很 高 的 
要 求 。 目前 的 SKA 探 路 者 望远镜 由 于 受到 灵敏 度 和 分 辩 率 的 限制 ， 在 宇宙 黎明 和 再 电离 探 
测 方面 主要 致力 于 统计 测量 ， 而 未 来 的 SKA 则 同时 具备 统计 测量 和 成 像 的 能 力 。 然 而 随 着 
SKA 灵敏 度 和 分 辩 率 的 提高 ， 为 满足 宇宙 黎明 和 再 电离 信号 探测 需求 所 需 采 集 的 数据 量 的 


增 大 ， 大 视 场 、 带 宽 、 高 分 辨 率 以 及 多 波束 技术 ， 对 硬件 和 软件 都 提出 了 更 高 和 更 有 具 挑战 性 
的 要 求 。 


在 低频 射电 波段 进行 探测 所 面临 的 主要 问题 之 一 是 干扰 去 除 问题 。 随 着 科技 的 发 展 ， 射 
电 王 扰 是 目前 低频 射电 探测 的 主要 困难 之 一 ， 射 电 宁 静 区 域 非 常 稀少 。 因 此 ， 射 频 干 扰 信 号 


的 识别 和 去 除 在 低频 射电 数据 处 理 中 非常 关键 。 低 频 射 电 干 扰 分 为 不 同类 型 ， 如 出 现在 某 个 
. 特定 频率 的 强 干 扰 ， 一 般 来 源 是 电视 基站 、 卫 星 信号 、 飞 机 、 雷 达 等 ， 出 现在 某 个 特定 时 间 
A 的 频率 覆盖 范围 较 宽 的 射电 干扰 ， 如 闪电 、 高 压 电 缆 等 。 对 于 这 些 射电 干扰 并 没有 统一 的 方 
CD 法 来 进行 扣除 ， 尤 其 在 面临 大 量 数 据 处 理 时 ， 需 要 更 为 有 效 的 自动 识别 和 去 除 干 扰 的 方法 。 
E 目前 能 够 对 观测 数据 进行 自动 干扰 信号 识别 和 去 除 的 方法 包括 对 REI 自动 化 处 理 的 软件 包 
c AOFLAGGER"", FLAGCAL""^&, AOFLAGGER 是 LOFAR 去 除 RFI 默认 的 检测 管道 ， 


是 在 时 频 空间 对 信号 进行 拟 合 ， 准 确 率 差 不 多 与 人 眼 效果 相当 。GMRT 主要 用 FLAGCAL 
进行 射电 干涉 测量 数据 的 RFI 标记 和 校准 ， 尽 管 在 本 质 上 与 其 他 算法 相似 ， 但 FLAGCAL 
拥有 强大 的 估算 系统 ， 并 且 可 以 做 到 高 性 能 计算 和 并 行 计算 。 
数据 校准 也 是 数据 处 理 中 的 重要 步骤 之 一 。 在 对 干涉 阵列 进行 数据 处 理 的 过 程 中 ， 只 有 
对 观测 数据 进行 校准 才能 获得 有 科学 价值 的 观测 结果 。 校 准 是 为 了 对 仪器 效应 以 及 电离 层 
影响 进行 去 除 ， 从 而 对 观测 的 射电 源 进行 流量 和 位 置 的 改正 。 在 进行 大 视 场 成 像 时 ， 源 于 射 
电源 的 信号 会 随 方向 、 时 间 和 频率 发 生变 化 ， 尤 其 对 于 大 视 场 成 像 而 言 ， 对 于 不 同 的 观测 方 
向 和 时 间 进 行 校准 非常 关键 。 数 据 校 准 过 程 分 为 幅度 校准 和 相位 校准 ， 即 拟 合 一 个 复数 的 增 
益 因 子 对 每 一 个 接收 单元 进行 改正 ， 对 观测 视 场 内 的 射电 源 进行 位 置 和 亮度 的 校准 。 在 校准 
过 程 中 ， 可 以 使 用 天 线 的 波束 模型 、 射 电源 天 空 模型 或 者 特定 的 校准 源 进行 数据 校准 。 在 实 
= 际 观 测 中 ， 精 准 的 天 线 波束 模型 的 建立 和 天 空 模型 的 建立 都 存在 困难 ， 因 此 校准 过 程 的 误差 
也 会 传递 到 最 终 的 宇宙 黎明 和 再 电离 信号 提取 的 过 程 中 。 
CLEAN 是 目前 普遍 使 用 的 成 像 方法 之 一 。 基 于 目前 低频 望远镜 阵列 数据 处 理 的 需求 ， 
目前 开发 出 的 成 像 软 件 包 包括 WSCLEAN (w-Stacking Clean) 等 ， 这 些 软件 包 都 具备 大 视 
场 成 像 中 w 项 的 改正 功能 。WSCLEAN” 是 用 C++ 语言 编写 的 ， 使 用 C++ 11 标准 库 的 
线程 模块 可 以 对 多 个 步骤 进行 多 线程 处 理 : 读 和 写 、 不 同 w 项 之 间 的 网 格 化 、 执行 FFTs 
以 及 Hógbom Clean" 4%. WSCLEAN 当前 使 用 Kaiser-Bessel (KB) 窗口 函数 ， 该 函数 易 
于 计算 且 快 速 ， 非 常 类 似 于 PSWF (the prolate spheroidal wave function). 


T 
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在 宇宙 再 电离 探测 中 ， 目 前 存在 的 数据 校准 和 功率 谱 计 算 软 件 主要 包括 the fast 
holographic deconvolution (FHD), eppsilon, the cosmological HI power spectrum estima- 
tor (CHIPS) “. FHD/eppsilon 功率 谱 分 析 ” 在 研究 中 很 普遍 ， 包 括 MWA Phase I”, 
MWA Phasell™, PAPER” 以 及 在 计划 中 的 HERA 的 数据 处 理 ， 也 可 用 于 仿真 数据 ”。 
FHD/eppsilon 结合 对 光谱 的 精确 校准 ， 以 及 端 到 端 数据 匹配 的 模型 仿真 和 误差 传递 ， 非 常 
适合 宇宙 再 电离 信号 的 测量 、 数 据 处 理 和 研究 。CHIPSW 是 MWA EoR 使 用 的 信号 处 理 管 
线 之 一 ， 用 于 获取 校准 后 的 数据 并 将 其 处 理 为 输出 功率 谱 。 全 部 热 噪声 加 上 残余 前 景 信 号 
的 道 协 方差 加 权 是 提取 信和 号 统计 信息 的 最 佳 方法 ， 但 这 存在 一 定 程度 的 困难 ” 。 残 差 前 景 
模型 不 恰当 ， 以 及 无 法 完全 独立 地 了 解 其 内 部 协 方 差 和 与 信号 的 协 方差 ， 容 易 导 致 信号 丢 
AU. pt, CHIPS 在 整个 工作 中 主要 被 用 作 逆 方差 估计 量 ， 其 中 基线 权重 用 于 采样 。 

测量 高 红 移 HI 信号 所 处 的 波段 是 低频 波段 50 ~ 200 MHz， 比 起 其 他 较 高 的 观测 频率 ， 
这 一 波段 存在 的 人 为 干扰 较 多 ， 去 除 和 识别 起 来 比较 困难 。 由 于 前 景 射电 源 暴 律 谱 的 特性 ， 


辐射 在 低频 波段 更 强 ， 也 造成 了 更 强 的 前 景 污染 。 并 且 ， 目 前 低频 观测 图 像 的 分 辨 率 和 图 像 
7 质量 低 于 其 他 波段 ， 造 成 了 前 景 射 电源 精确 识别 的 限制 和 建立 精确 前 景 模型 的 限制 ， 而 数据 
> BEE LAE T-ALL ROT TFS Ak BL SA IE AN 4e 
po 准确 ， 误 差 会 造成 数据 校准 的 困难 。 同 时 ， 在 低频 波段 的 电离 层 影响 会 引起 射电 源 位 置 和 亮 
Co 度 的 变化 ， 在 分 辩 率 受 限 的 情况 下 ， 校 准 也 存在 困难 。 


在 宇宙 再 电离 功率 谱 探测 的 过 程 中 ， 最 理想 的 处 理 方 式 是 使 用 原始 的 观测 数据 直接 构 
造 功率 谱 ， 避 免 格 点 化 问题 引起 的 误差 。 然 而 低频 射电 阵列 原始 观测 的 数据 量 巨大 。 另 外 为 
达到 提取 微弱 信号 的 需求 ， 需 要 进行 长 时 间 的 数据 积累 。 并 且 要 完成 对 原始 数据 的 校准 ， 对 
计算 和 存储 有 很 大 需求 。 因 此 在 宇宙 再 电离 探测 的 数据 处 理 中 面临 着 海量 数据 处 理 的 问题 。 

未 来 SKA 将 采用 深 场 定点 观测 模式 和 巡天 模式 来 完成 对 宇宙 黎明 和 再 电离 时 期 深 场 成 
像 和 统计 测量 工作 。 大 天 区 的 成 像 涉及 到 w 项 改正 以 及 图 像 拼接 的 问题 ， 同 时 要 保证 提取 
言 号 所 需 的 高 动态 范围 。 因 此 大 视 场 成 像 也 是 数据 处 理 中 的 难点 之 一 。 


3 宇宙 黎明 和 再 电离 时 期 探测 进展 


在 过 去 10 年 间 国 际 上 已 经 有 许多 低频 探测 设备 建成 ， 人 们 利用 这 些 探测 设备 在 宇 
明和 再 电离 探测 领域 取得 了 研究 进展 。 
3.1 ”全 天 总 功率 测量 

宇宙 黎明 和 再 电离 时 期 的 全 天 总 功率 测量 通常 是 由 单 天 线 完 成 ， 测 量 方法 类 似 于 宇 
宙 微 波 背景 辐射 测量 。 目前 国际 上 的 相关 实验 包括 EDGES, BIGHORN, SCI-HI, LEDA, 
SARAS, DARE, PRIZM, ASSASSIN, DAPPER 和 DSL 等 。 

EDGES 实验 位 于 澳大利亚 西部 ， 则 在 测量 再 电离 信号 的 全 天 总 功率 ， 它 采用 了 一 种 
特定 的 方法 研究 预期 的 频谱 差异 ， 以 克服 前 景 和 信号 分 离 上 的 困难 ， 这 种 方法 是 限制 实验 范 
围 以 测试 光谱 中 的 不 连续 特征 ， 因 为 这 些 特征 必定 是 由 于 红 移 21 cm 谱 线 亮 温 度 的 快速 跃 
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迁 而 不 是 较 平滑 的 前 景 。2010 年 ，Bownman 和 Rogers" 率先 对 宇宙 再 电离 时 标 进行 了 约 
WR: 2018 年 ，Bownman 和 Rogers 发 表 了 其 在 低频 波段 全 天 总 功率 探测 的 实验 结果 ， 经 
过 数据 校准 ， 前 景 去 除 等 数据 处 理 过 程 ， 在 78 MHz 附近 检测 到 一 个 吸收 谷 ， 提 取出 疑似 
宇宙 黎明 时 期 的 HI 21 cm 辐射 的 吸收 信号 ;2017 年 ，Monsalve 等 人 "通过 高 频段 天 线 数 
据 ， 对 字 宙 再 电离 模型 进行 了 约束 。BIGHORN 特点 是 便捷 、 简 单 、 低 功 耗 ， 可 自由 选择 
试验 地 点 ，2018 Æ, Choudhury 和 Datta’ ” 提出 由 artificial neural networks (ANN) 从 明亮 
的 银河 系 前 景 中 探测 21 cm 全 天 信号 (以 “tanh” 模 型 的 形式 表示 )， 获 得 初步 结果 ， 发 现 即 
使 存在 由 三 阶 或 更 高 阶 对 数 多 项 式 表示 的 明亮 前 景 ， 也 可 以 准确 确定 全 天 21 cm 信号 参数 。 
SCLHIS 用 4.4 h 的 数据 ， 通 过 多 项 式 扣除 前 景 ， 在 60~88 MHz 范围 内 给 出 了 1 K 均 方 根 
IAWR. LEDA 用 19 min 的 有 效 测量 数据 ， 对 50~100 MHz 之 间 的 吸收 模 进 行 了 限 
制 。SARAS 位 于 印度 ， 是 一 个 单 天 线 相关 频谱 仪 ，Patra 45 AP 提供 了 110~175 MHz i 
段 射电 频谱 的 绝对 天 空 亮 度 和 谱 指 数 ， 并 对 150 MHz 全 天 图 提供 了 更 精确 的 校准 。SARAS 
27" seus BEER T Cohen 等 人 "模拟 的 几 个 再 电离 时 期 全 天 21 cm 信号 模型 。 地 面 的 这 些 实 
验 ， 虽 然 都 选择 了 人 烟 稀少 、 电 磁 环 境 优 良 的 测量 地 点 ， 但 在 一 定 程度 上 仍 受到 人 为 电磁 
干扰 的 影响 ， 并 且 也 会 受到 地 球 电离 层 吸 收 、 折 射 的 影响 ， 如 果 能 在 太空 进行 观测 ， 特 别 
是 利用 月 球 侦 挡 来 自 地 球 的 电磁 干扰 ， 将 大 大 提升 实验 的 灵敏 度 和 可 靠 性 ， 美 国 为 此 提出 
T DARE 实验 ，DARE 太空 计划 打算 从 月 球 远 端 上 方 的 轨道 测量 全 天 21 cm 信号 ”， 通 过 
它 在 黑暗 时 期 末 和 宇宙 黎明 期 间 对 银河 系 星际 介质 的 影响 ， 来 追踪 第 一 代 发 光 天 体 的 形成 ， 
DARE 绕 月 球 飞行 的 任务 为 期 3 年 ， 并 从 月 球 上 方 获取 数据 ， 这 是 太阳 系 中 唯一 经 证 实 不 
受 射频 干扰 和 任何 重要 电离 层 影 响 的 位 置 。 PRIZSM 是 测量 全 天 平均 天 空 亮 度 的 新 实验 ， 其 
中 包括 由 第 一 代 恒 星 形成 引起 的 预期 的 红 移 21 cm HI 吸收 特征 ，PRIZM 由 两 个 双 极 化 天 
线 组 成 ， 工 作 的 中 心 频率 分 别 为 70 MHz 和 100 MHz， 该 实验 位 于 亚 南极 的 马里 恩 岛 上 ， 马 
里 恩 岛 拥有 异常 干净 的 观测 环境 ，FM 频段 内 基本 上 没有 可 见 的 污染 。ASSASSIN BERRA 
宇宙 黎明 到 再 电离 时 期 结束 这 一 阶段 的 红 移 21 cm 信号 。DAPPER 主要 是 研究 早期 宇宙 的 
黑暗 时 期 ， 是 一 项 绕 月 轨道 实验 ， 这 艘 飞船 将 使 用 专门 的 接收 器 和 高 频 天 线 来 监听 来 自 早 
期 宇宙 的 微弱 信号 ， 要 做 到 这 一 点 需要 一 个 非常 安静 的 环境 ， 而 月 球 的 背面 是 一 个 理想 的 
地 方 。2014 年 ， 中 欧 联 合 团队 提出 了 DSL 计划 ， 对 超 长 波 天 文 观测 阵列 方案 开展 了 系统 研 
= 究 ， 主 要 是 月 球 轨道 上 的 线性 阵列 ， 同 时 对 天 空 进行 干涉 测量 ， 可 以 合成 绘制 整个 天 空 的 高 
w 分 辨 率 图 像 以 及 全 天 射电 频谱 ， 研 究 低 频 射 电源 、 银 河 系 星际 介质 等 。 全 天 总 功率 测量 实验 
详细 的 信息 及 进展 情况 如 表 1 所 示 。 
目前 的 SKA 探 路 者 阵列 ， 一 般 采 取 定 点 观测 和 巡天 观测 两 种 模式 。 相 对 于 LOFAR， 
MWA 等 采用 定点 观测 和 巡天 观测 模式 的 阵列 ， 中 国 的 21CMA 阵列 的 天 线 是 单一 指向 天 
线 ， 视 场 是 北 天 极 区 域 ， 优 势 是 可 以 全 天 指向 北 天 极 区 域 实现 长 时 间 的 数据 积累 ， 缺 点 是 受 
到 了 观测 天 区 的 限制 ， 在 定 标 源 选取 上 和 展 源 的 去 除 上 会 遇 到 困难 。 

针对 全 天 21 cm 信号 ，2017 Æ Jordan 等 人 ”利用 星系 和 暗 举 的 丰 度 匹配 建立 了 一 个 新 
模型 ， 把 21 cm 信号 和 高 红 移 星系 光度 函数 联系 起 来 。 该 模型 已 针对 高 红 移 星系 光度 函数 的 
测量 进行 了 校准 ， 并 进一步 调整 以 匹配 宇宙 微波 背景 汤姆 逊 散射 的 光 深 。 模 型 在 110 MHz 
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表 1 全 天 总 功率 测量 实验 


频率 f 


i AUR 参考 文献 
50-100 6.5~14.8 BH 
XG 电离 时 标 、 模 型 进行 了 约束 ， 提 出 疑似 宇 
EDGES 和 和 iei i a, 时 标 、 模 型 进行 了 约束 ， 提 出 疑似 
(HES. l | 


100—200 13.2~27.4 


ANN 从 明亮 的 银河 系 前 景 中 探测 21 全 天 信 
BIGHORN 20~160 7.9~70 4. ee ps UAR TEL 2T GEB ele fr] 
号 (以 “tanh” 模 型 的 形式 表示 )。 


用 4.4 h 的 数据 ， 通 过 多 项 式 扣 除 前 景 ， 在 60~88 

SCLHI 40-130 9.9~34.5 g ne pu dat En bans maj 
MHz 范围 内 给 出 了 1 K 均 方 根 残 差 的 限 值 。 

用 19 min 的 有 效 测量 数据 ， 对 50~100 MHz 之 间 

LEDA 50~100 — 13.2~27.4 UAM : | 
的 吸收 模 进 行 了 限制 。 Ira, T3 
SARAS 提供 了 110—175 MHz 波段 的 绝对 天 空 亮 

SARAS ”87.5~175 7.1~15.2 ” 度 和 谱 指数 ，SARAS 2 排除 了 Cohen 等 人 模拟 的 E] 


几 个 再 电离 时 期 全 天 21 cm 信号 模型 。 


IT 


x DARE 40~120 ”10.8~34.5 计划 从 月 球 远 端 上 方 的 轨道 测量 全 天 21 cm 信号 。 (22, 23] 
y 70 和 描述 了 其 最 初 的 设计 和 配置 ， 在 极 宁静 的 射电 场 演 
PRESE 400 示 了 全 天 21 cm 探测 设备 的 运行 情况 。 pal 
采用 了 紧密 排列 的 天 线 阵列 而 不 是 使 用 传统 方法 使 


= 用 单个 天 线 来 测量 。 用 仿真 演示 了 最 新 开发 的 广 
e ASSASSIN n2 127.4 
© SSASSIN 50~200 6I~274 法， 并 表明 对 于 理想 化 的 仪器 ， 理 论 上 可 以 从 紧密 pa 


间隔 的 偶 极 子 阵列 的 可 见 性 中 提取 21 cm 的 信号 。 
BABBER- dees € 将 测量 光谱 的 幅度 ， 使 POS 262 eon) 标 四 
准 宇宙 学 模型 与 其 他 标准 宇宙 模型 所 需 的 水 平 。 

这 个 项 目的 主要 目标 是 使 用 卫星 星座 作为 干涉 阵列 

来 绘制 30 MHz 以 下 的 天 空 谱 ， 另 一 个 目标 是 在 感 
DSL 30~100 ”13.2~46 ”兴趣 的 频率 范围 内 进行 高 精度 的 全 天 频谱 测量 ， 以 Ez 
= 研究 黑暗 时 期 和 宇宙 黎明 ， 因 为 不 受 地 面 观 测 中 的 
gm 电离 层 干 扰 和 射频 干扰 。 


处 吸收 达到 峰值 。 若 不 引入 新 物理 模型 ， 无 法 解释 EDGES 测量 到 的 疑似 信号 。 如 图 3 所 
示 。 相 较 之 下 ，EDGES 得 到 的 全 天 信号 ”的 确 存 在 一 个 吸收 谷 ， 这 符合 再 电离 时 期 前 的 
冷 IGM 模型 ， 但 同时 也 有 一 些 新 的 特征 ， 如 : 信号 的 吸收 峰值 达到 了 约 0.5 K， 这 远 远 超 
出 了 此 前 的 预测 ， 峰 值 频率 在 78 MHz， 亦 即 更 高 的 红 移 。2018 E, Jordan 等 人 号 对 此 差 
异 给 出 解释 : (1) 吸收 峰值 的 提前 可 能 源 自 于 高 红 移 矮星 系 恒星 形成 效率 的 提高 ， 从 而 导致 
Wouthuysen-Field 效应 和 X 射线 加 热 的 提前 ; (2) 建立 一 个 热 历史 演化 的 参数 模型 ， 提 供 
IGM 新 的 冷却 途径 ， 来 拟 合 EDGES 的 深 吸 收 谷 ; (3) 观测 到 深 吸 收 谷 也 说 明 ， 除 了 宇宙 背 
景 辐射 以 外 ， 可 能 还 存在 新 的 射电 背景 如 恒星 形成 星系 提供 的 射电 源 ， 利 用 已 有 恒星 形成 速 
率 与 射电 光度 的 经 验 关 系 ， 可 以 对 EDGES 观测 做 出 拟 合 。 
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AT, /mK 


一 - Bowman 等 人 2018 
-一 Jordan 等 人 2017 


0 50 100 150 200 
v/MHz 


He 图 中 红线 为 Jordan 等 人 中 针对 全 天 21 cm 信号 建立 的 新 模型 ， 可 以 为 21 cm 信号 和 高 红 移 星系 光度 函数 


建立 联系 。 蓝 色 表 示 基 于 Bownman 和 Rogers 的 观测 结果 的 拟 合 结果 (Flatten Gaussian), TE 78 MHz 


附近 检测 到 一 个 吸收 谷 ( 浅 灰 


3.2 ”功率 谱 测量 


色 虚 线 )， 疑 似 宇宙 黎明 时 期 信号 。 


3 ”针对 21 cm 信号 建立 的 模型 


目前 国际 上 的 SKA 探 路 者 低频 射电 阵列 ， 如 21CMA, MWA, LOFAR, GMRT, PAPER, 


HERA 等 ， 都 将 宇宙 黎明 和 再 电离 时 期 的 功率 谱 测量 作为 首要 科学 目标 之 一 。 宇 宙 黎 明和 


再 电离 时 期 的 功率 


普 测 量 分 为 两 个 方向 : 通过 去 除 射 电 前 景 污 染 ， 构 建 宇 宙 再 电离 背景 信号 


的 功率 谱 ;， 直接 构建 功率 谱 ， 通 过 在 功率 谱 空 间 研究 宇宙 再 电离 信号 的 探测 窗口 ， 避 免 仪 器 


效应 和 前 景 污染 的 


区 域 ， 直 接 提取 宇宙 再 电离 时 期 信号 。 近 年 来 ， 基 于 这 些 低 频 射电 阵列 的 


观测 数据 ， 在 统计 测量 方面 已 经 进行 了 大 量 的 研究 工作 ， 如 表 2 所 示 。 


表 2 ”功率 谱 统计 测量 实验 ( 表 中 仅 罗 列 出 已 有 测量 结果 的 实验 ) 


实验 频率 f 红 移 > RE k 测量 结果 参考 
名 称 /MHz /(h-Mpc^!) 03K? 文献 
175.3 7.1 0.27 A? <2.7 x10* BZ 
177.5 7.0 0.20 A? <3.9 x10° 

MWA 7 " am 
189.3 6.5 0.59 A? «2.39 x10 Ba 
189.3 6.5 0.14 (43 mK)? = 1.8 x 10° [a] 
127.9 10.1 e 0.053 A? < 79.6? 

LOFAR " à 
22140.6 z 9.1 0.075 ^? < 73 Ba 
GMRT 147.9 8.6 0.5 A? « 248? Bo 
151.1 8.4 0.15~0.5 A? < 22.4? BI 

PAPER " 2 
151.5 8.37 0.37 A? < 200 [EU] 
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为 进行 前 景 射电 源 属性 研究 ， 建 立 精 确 的 天 空 模型 ， 目 前 基于 低频 射电 阵列 的 观测 数 
据 ， 已 经 形成 的 射电 源 星 表 包 括 : 21CMA Rx". MWA GLEAM ER 等。 

在 功率 谱 测 量 方面 ，Procopio 等 人 四 基于 MWA 的 GLEAM 巡天 ， 为 提高 分 辨 率 ， 结 
A GMRT TGSS 巡天 150 MHz 的 数据 ， 为 致密 的 前 景 源 建立 了 一 个 模型 ， 这 种 混合 模型 可 
以 对 EoR1 (MWA EOR 实验 的 三 个 目标 场 之 一 ，(RA, DEC) = (4 h, —30.0°)) 的 数据 进行 重 
新 处 理 ， 最 大 可 提高 8% 的 校准 精度 。 利 用 该 模型 ，Beardsley AE k = 0.27 h-Mpc7}, 
:二 7.1 得 到 功率 谱 A? «27 K^; Barry ATE k = 0.2 h-Mpc-!, z = 7.0188] A? <3.9 K?; 
Li & A" 在 大 = 0.59 h-Mpc!, z = 6.5 得 到 A2 <2.39 K?; Trott & A" RIE MWA EoR 项 
目的 数据 ， 通 过 实时 系统 的 校准 ， 以 及 CHIPS 管线 估算 功率 谱 ， Æ k = 0.07 ~ 3 h-Mpc-!, 
2=65~ 8.7 计算 出 了 平均 功率 谱 ，Byrne” 提出 了 南部 天 空 的 全 极 化 漫 射 图 ， 并 讨论 
了 其 用 于 改进 EoR 校准 和 前 景 去 除 的 应 用 。Patil 5& AP" YE k © 0.053 h-Mpc-!, z = 10.1 
得 到 A? < 79.6? x 10-3K?; Mertens 4& AP" 基于 LOFAR 测 得 的 141 h 的 数据 ， 获 得 了 在 
z 9.1 的 红 移 下 21 cm 信号 功率 谱 新 的 测量 上 限 ， 这 项 研究 在 频谱 平滑 增益 校准 ， 高 
斯 过 程 回 归 (GPR) 前 景 缓 解 和 最 佳 加 权 功 率 谱 推断 等 方面 都 进行 了 重大 改进 ，Bradley 
等 人 ™ 用 21CMMC 对 3D 宇宙 21 cm 信号 进行 正 向 建 模 ， 探 索 了 与 观测 极限 不 同 的 天 体 
物理 参数 空间 区 域 。PAPER 项 目 组 Ali ATE k = 0.15 ~ 0.5hMpc-1, z = 8.4 得 到 
A? < 22.42 x 10-3K?; Cheng 等 人 ”通过 对 其 数据 进行 可 靠 的 重新 分 析 ， 改 进 了 之 前 的 测 
量 结果 ，Kolopanis 等 人 ”由 PAPER 的 数据 ， 给 出 了 再 电离 时 期 21 cm 功率 谱 的 限制 ， 该 
工作 对 数据 处 理 管 线 中 的 基于 延迟 的 前 景 滤波 以 及 基于 经 验 的 协 方差 估计 等 步 又 进行 了 改 
进 ， 并 且 解 决 了 PAPER 先前 的 结果 中 的 信号 损失 问题 ，Jacobs 等 人 ”基于 PAPER 数据 ， 
发 表 了 145 MHz 频率 的 功率 谱 测 量 结果 ， 该 测量 有 望 涵盖 宇宙 电离 频段 内 的 大 部 分 信号 
源 ， 在 建立 前 景 信号 模型 方面 取得 了 进展 。GMRT 目前 也 在 进行 高 红 移 HI 功率 谱 的 测量 工 
作 ，Paciga 等 人 "用 奇异 值 分 解 (SVD) 的 方法 来 去 除 前 景 ， 对 前 景 结 构 的 假设 条 件 引 入 较 
少 。 通 过 这 种 方法 ， 他 们 还 首次 量化 了 由 前 景 滤波 器 引起 的 信号 损耗 ， 并 提出 了 针对 该 损 
耗 调 整 的 功率 谱 测量 ， 从 而 得 到 了 在 有 = 0.5 h:Mpc-1，z = 8.6 时 的 2c 上 限 的 修正 测量 值 
A? < 2482 x 10-3K2。 

我 们 对 低频 射电 阵列 在 宇宙 黎明 时 期 功率 谱 研 究 方面 的 最 新 探测 结果 与 理论 预测 信和 号 
强度 进行 了 比较 ， 结 果 如 图 4 所 示 。Matias A" 根据 理论 模型 得 出 21 cm 信号 的 峰值 在 
0.01 K AA, WA 4 中 黑色 虚线 所 示 。 

3.3 ”成 像 研究 

未 来 的 SKA 低频 阵列 不 仅 能 对 宇宙 再 电离 信号 进行 功率 谱 测量 ， 还 将 具备 对 宇宙 再 电 
离 信号 进行 成 像 研究 的 能 力 ， 完 成 低频 射电 干涉 大 视 场 、 高 动态 、 多 波束 的 成 像 处 理 软件 ， 
将 有 助 于 揭 开 宇宙 再 电离 的 神秘 面纱 。 与 功率 谱 测量 相 比 ， 对 再 电离 时 期 的 大 尺度 电离 氢 区 
进行 直接 成 像 观测 ， 可 以 更 为 直接 地 提取 电离 区 域 的 大 小 和 特征 等 物理 信息 。 对 5 个 20 F 
方 度 的 天 区 进行 深度 成 像 观 测 也 是 未 来 SKA 低频 的 首要 科学 目标 之 一 。 

SKA 是 由 全 球 超过 10 个 国家 合资 建造 的 世界 最 大 综合 孔径 射电 望远镜 阵列 ， 台 址 位 于 
澳大利亚 、 南 非 及 南部 非洲 8 个 国家 的 射电 宁静 区 域 ， 由 分 布 在 3000 km 范围 内 的 约 2 500 
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---- AT=0.01 K 
e MWA 
e LOFAR 
e GMRT 
€ PAPER 


AT /mK 


Ha 红色 是 MWA， 绿 色 是 LOFAR， 蓝 色 是 GMRT， 紫 色 是 PAPER. SEE 20. Matias 等 人 四 根据 
理论 模型 得 出 21 cm 信号 的 峰值 在 0.01 K 左右 ， 如 黑色 虚线 所 示 。 


4 ”功率 谱 测 量 结果 汇总 (与 表 2 对 应 一 致 ) 


面 15 m 口径 碟 形 高 频 天 线 、250 km 范围 内 的 250 个 直径 约 60 m 的 致密 孔径 阵列 (中 频 ) 以 
及 130 万 个 对 数 周 期 天 线 组 成 的 稀 朴 孔径 阵列 (低频 ) 组 成 ， 其 等 效 接收 面积 达 平 方 公里 级 ， 
频率 覆盖 范围 为 50 MHz ~ 20 GHz. SKA 分 为 建设 准备 阶段 (2012 一 2021 年 )、 建 设 第 一 阶 
Be (SKA1, 2021—2028 年 ) 和 建设 第 二 阶段 (SKA2，2028 年 后 )。SKA 的 首要 科学 目标 之 

就 是 利用 HI 21 cm 信号 进行 EoR 成 像 观测 ， 以 及 进行 CD/EoR 功率 谱 测 量 。 与 目前 正在 
运行 的 SKA 探 路 者 项 目 相 比 ，SKAI 低频 阵列 是 唯一 可 实现 EoR 成 像 观 测 的 实验 装置 ， 具 
重大 的 历史 意义 。SKA 低频 阵列 的 大 面积 巡天 还 可 以 统计 测量 红 移 6 ~ 27 的 21 cm 信和 号 
功率 谱 ， 给 出 了 气体 的 三 维 分 布 ， 其 峰值 信 噪 比 约 为 100， 波 数 太 = (0.02 ~ 0.1) -Mpc ^ !. 
未 来 ，SKA2 将 大 大 提高 角 分 辨 率 和 灵敏 度 ， 在 更 大 的 尺度 范围 内 对 CD/EoR 进行 更 精确 
的 测量 。 由 于 宇宙 21 cm 信和 号 非常 微弱 ， 因 此 需要 进行 长 时 间 的 观测 来 降低 噪音 。 利 用 
SKA 探 路 者 低频 射电 阵列 进行 某 些 特定 天 区 的 深 场 观测 ， 可 以 更 好 地 了 解 前 景 的 性 质 、 仪 
器 的 影响 ， 改 善 数 据 处 理 的 方法 等 ， 同 时 也 为 准确 地 在 二 维 功率 谱 空间 限制 宇宙 再 电离 探测 
的 窗口 提供 了 可 能 性 。 这 些 对 于 在 不 和 久 的 将 来 使 用 SKA 低频 阵列 探测 宇宙 再 电离 信号 是 至 
关 重 要 的 。SKA 预期 在 100 ~ 200 MHz 获得 第 一 幅 宇宙 再 电离 的 直接 图 像 ， 揭 开 宇 宙 再 电 
离 的 奥秘 ， 为 了 这 一 极 富 挑 战 性 的 观测 ，Zheng 等 人 "对 现 有 的 目录 和 频率 范围 内 的 观测 
进行 了 预选 ， 选 出 适 于 宇宙 黎明 和 再 电离 深度 成 像 的 “最 安静 ”和 “最 干净 ”的 候选 天 区 。 
SKA1 的 基本 参数 概括 在 表 3 中 。 目 前 与 SKA 轿 盖 频段 相似 的 射电 望远镜 的 观测 灵敏 度 对 


= 


比如 图 5 所 示 ， 图 来 自 国际 SKA 组 织 (SKA Organization, SKAO™). 


V https: //www.skatelescope.org/key-documents/SKA Phase 1 Construction Proposal 
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#3 SKA 主要 技术 指标 


参数 SKA1 - low SKA1 - mid 
频率 范围 /GHz 0.05 ~ 0.35 0.35 ~ 20 
基准 频率 /GHz 0.11 1.67 
灵敏 度 (Aca /Tsys)/(m K7!) (基准 频率 ) 550 1500 (包括 MeerKAT) 
视 场 /(*)? (基准 频率 ) 14 0.33 
最 高 分 辨 率 /(”) (基准 频率 ) 7 0.25 
单元 阵列 主 频 范围 /GHz 0.1 ~ 0.2 0.35 ~ 20 
单元 阵列 尺寸 /m 35 15 
单元 阵列 数目 512 (规划 )，476 (实施 ) ”133 十 64 (MeerKAT)，130 (实施 ) 
每 单元 阵列 天 线 数 256 1 


最 长 基线 /km 


65 (规划 )，40 (实施 ) 


154 (规划 )，120 (实施 ) 


[T TTTTTT T T TTTTTIT 


T TT 


SKA1-LOW 
1 000 : 
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Cos rrrul 


频率 /GHz 


5 ， 目前 与 SKA 覆盖 频段 相似 的 射电 望远镜 


的 观测 灵敏 度 (来 自 SKAO) 


高 动态 大 视 场 成 像 是 未 来 SKA 低频 成 像 数 据 处 理 中 面临 的 主要 问题 之 一 。 为 满足 未 


K SKA 大 视 场 成 像 的 运算 需求 ， 据 估计 ， 计 算 机 硬 伯 
频 预算 中 占有 


设施 建设 所 产生 的 费用 将 在 SKA 低 


半 的 比例 。 在 处 理 未 来 SKA 大 视 场 成 像 的 过 程 中 ， 在 校准 、 高 动态 、 频 率 
带宽 、 天 线 波束 指向 、 主 斩 边 缘 、 电 离 层 改正 以 及 定 标 源 选取 等 问题 上 ， 都 存在 困难 与 挑 
战 。 大 视 场 范围 内 包含 了 大 量 的 射电 源 以 及 来 自 仪 器 、 大 气 层 的 影响 ， 都 随 方向 、 


时 间 、 


频率 的 变化 而 变化 ， 因 此 对 不 同 观测 方向 进行 有 效 的 校准 显得 极为 重要 。 目 前 ， 基 于 和 剥离 
(peeling)” 概念 的 校准 方法 是 唯一 有 效 的 ， 但 仍 需 继续 开发 从 而 降低 成 本 并 提高 计算 精度 ; 
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高 质量 的 图 像 在 宇宙 黎明 和 再 电离 时 期 的 数据 分 析 和 处 理 过 程 中 起 到 很 大 的 作用 : 建立 模 
型 、 校 准 数据 、 提 取 再 电离 信号 和 去 除 前 景 源 。 对 于 SKA 这 类 大 视 场 望远镜 ， 在 成 像 过 程 
中 需要 根据 天 线 的 不 同 指向 进行 相位 的 改正 ;给 定 的 带宽 情况 下 ， 积 分 时 间 决 定 了 图 像 信 噪 
比 ， 积 分 时 间 的 选择 决定 了 亮 的 定 标 源 在 保证 成 像 不 被 污染 时 它们 之 间 的 距离 。 


A 总 结 与 展望 


宇宙 黎明 和 再 电离 时 期 的 探测 是 通过 测量 HI 超 精细 能 级 跃迁 产生 的 21 cm 信号 来 研究 
第 一 代 发 光 天 体 的 产生 和 形成 过 程 。 由 于 宇宙 黎明 和 再 电离 时 期 的 HI 21 cm 信和 号 非常 微弱 ， 
因此 在 前 景 去 除 以 及 数据 处 理 的 各 个 环节 中 都 存在 着 困难 和 挑战 。 宇 宙 黎 明和 再 电离 时 期 
探测 主要 分 三 个 方向 : 全 天 总 功率 测量 、 功 率 谱 测量 以 及 直接 成 像 研究 。 全 天 总 功率 测量 由 
单 天 线 完 成 ， 功 率 谱 测量 和 直接 成 像 研究 由 低频 射电 阵列 完成 。 低 频 射电 阵列 可 以 通过 漂 
移 扫描 模式 ， 进 行 浅 度 巡 天 ， 实 现 大 视 场 成 像 ， 在 50 ~ 200 MHz 低频 波段 ( 即 红 移 范围 6 
~ 27) 统计 获得 宇宙 黎明 和 再 电离 时 期 的 HI 功率 谱 。 宇 宙 黎 明和 再 电离 时 期 的 功率 谱 测量 
分 为 两 个 方向 : 通过 去 除 射 电 前 景 污染 ， 构 建 宇 宙 再 电离 背景 信号 的 功率 谱 ; 直接 构建 功率 
谱 ， 通 过 在 功率 谱 空 间 研究 宇宙 再 电离 信号 的 探测 窗口 ， 避 免 仪 器 效应 和 前 景 污染 的 区 域 ， 
直接 提取 宇宙 再 电离 时 期 信号 。 功 率 谱 测量 是 目前 许多 SKA 低频 探 路 者 阵列 的 首要 科学 目 
标 之 一 ， 由 于 受到 灵敏 度 和 分 辨 率 的 限制 ， 目 前 的 SKA 探 路 者 阵列 不 具备 直接 成 像 观测 的 
能 力 。 未 来 的 SKA 低频 阵列 将 进行 宇宙 再 电离 时 期 功率 谱 测 量 和 成 像 研究 ， 有 具备 对 宇宙 的 
再 电离 区 域 直 接 成 图 的 能 力 。 与 功率 谱 测量 相 比 ， 对 再 电离 时 期 的 大 尺度 电离 氨 区 进行 直接 
成 像 观 测 ， 可 以 更 为 直接 地 提取 电离 区 域 的 大 小 和 特征 。 本 文 主要 讨论 了 在 宇宙 黎明 和 再 电 
离 时 期 探测 中 关于 前 景 去 除 、 数 据 处 理 方法 等 相关 问题 ， 以 及 低频 射电 望远镜 在 全 天 总 功率 
测量 、 功 率 谱 测 量 方面 的 研究 进展 。 

低频 射电 阵列 在 灵敏 度 和 分 辨 率 上 的 提高 ， 以 及 随 着 软件 和 硬件 发 展 带 来 的 数据 处 理 
水 平 的 提高 ， 使 最 终 探 测 到 宇宙 黎明 和 再 电离 信号 的 可 能 性 也 不 断 增 加 。 宇 宙 黎 明和 再 电离 
时 期 探测 对 于 理解 宇宙 中 第 一 代 发 光 天 体 是 如 何 形成 的 和 何 时 形成 的 ， 电 离 氢 区 的 尺度 和 
成 团 性 等 特征 ， 宇 宙 再 电离 时 期 的 关键 物理 过 程 ， 以 及 宇宙 早期 结构 形成 乃至 整个 星系 的 形 
成 和 演化 历史 都 具有 重要 意义 。 对 宇宙 黎明 和 再 电离 信号 的 功率 谱 测量 和 成 像 研究 ， 将 揭示 
宇宙 从 黑暗 走向 光明 的 历史 ， 以 及 宇宙 结构 形成 从 线性 向 非 线性 的 演化 进程 。 
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Progress on the Exploration of the Cosmic Dawn and 


the Epoch of Reionization 


ZHAO Bi-xuan!, ZHENG Qian'?, ZHANG Ze-kang!?, 
GUO Quan!?, SHAN Huan-yuan!?? 


(1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 
2. Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210033, China; 3. Uni- 
versity of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: The cosmic dawn (CD) and epoch of reionization (EoR) are two of the most 
important missing pictures in the history of the Universe. The first luminous objects were 
formed during the CD/EoR period and thus marked the end of the dark ages. Detecting 
the redshifted 21cm hyperfine transition line of neutral hydrogen is regarded as the most 
promising and effective one to probe the CD/EoR. It includes three aspects: the total pow- 


er measurement of the CD/EoR, the statistical measurement of the CD/EoR signal and 
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perform imaging observations on the reionization area directly. Recently, a lot of related 
works have been done on both theories and observations, including i) theories and numerical 
simulations on the IGM evolution history during EoR; ii) data processing methods includ- 
ing radio frequency interference (RFI) removal, data calibration, techniques of foregrounds 
removal, etc. It makes great progress on the global detections and statistic measurements of 
the CD/EoR signals. As the largest radio telescope ever built, the Square Kilometer Array 
(SKA) will enable us a better understanding of the history and morphology of the CD/EoR. 
This work will provide a brief overview of the related research progress of CD/EoR detection 


in recent years. 


Key words: cosmology; the cosmic dawn and the epoch of reionization; radio interference 
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